A Runge-Kutta Method for Index 1 Stochastic Differential-Algebraic Equations with Scalar Noise
نویسندگان
چکیده
The paper deals with the numerical treatment of stochastic differential-algebraic equations of index one with a scalar driving Wiener process. Therefore, a particularly customized stochastic Runge-Kutta method is introduced. Order conditions for convergence with order 1.0 in the mean-square sense are calculated and coefficients for some schemes are presented. The proposed schemes are stiffly accurate and applicable to nonlinear stochastic differential-algebraic equations. As an advantage they do not require the calculation of any pseudo-inverses or projectors. Further, the mean-square stability of the proposed schemes is analyzed and simulation results are presented bringing out their good performance.
منابع مشابه
Projected Runge-kutta Methods for Differential Algebraic Equations of Index 3
In the present paper we introduce a new class of methods, Projected RungeKutta methods, for the solution of index 3 differential algebraic equations (DAEs) in Hessenberg form. The methods admit the integration of index 3 DAEs without any drift effects. This makes them particularly well suited for long term integration. Finally, implemented on the basis of the Radau5 code, the projected Runge-Ku...
متن کاملNumerical Solution of Stochastic Differential Equations with Constant Diffusion Coefficients
We present Runge-Kutta methods of high accuracy for stochastic differential equations with constant diffusion coefficients. We analyze L2 convergence of these methods and present convergence proofs. For scalar equations a second-order method is derived, and for systems a method of order one-and-one-half is derived. We further consider a variance reduction technique based on Hermite expansions f...
متن کاملLagrange-d’Alembert SPARK Integrators for Nonholonomic Lagrangian Systems
Lagrangian systems with ideal nonholonomic constraints can be expressed as implicit index 2 differential-algebraic equations (DAEs) and can be derived from the Lagrange-d’Alembert principle. We define a new nonholonomically constrained discrete Lagrange-d’Alembert principle based on a discrete Lagrange-d’Alembert principle for forced Lagrangian systems. Nonholonomic constraints are considered a...
متن کاملNumerical Methods for Second-Order Stochastic Differential Equations
We seek numerical methods for second-order stochastic differential equations that reproduce the stationary density accurately for all values of damping. A complete analysis is possible for scalar linear second-order equations (damped harmonic oscillators with additive noise), where the statistics are Gaussian and can be calculated exactly in the continuous-time and discrete-time cases. A matrix...
متن کاملEconomical Runge-Kutta Methods with Weak Second Order for Stochastic Differential Equations
In recent years there has been much interest in designing numerical methods for stochastic differential equations (SDEs) since they are used for the description of many real-life phenomena in different fields, including biology and physics, population dynamics, economics and finance. Because of the difficulty of obtaining exact solutions to such equations, numerical methods are needed. In this ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2010